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Accurate observations of surface ocean vector winds (OVW) with high spatial and 
temporal resolution are required for understanding and predicting tropical cyclones. As 
NASA’s QuikSCAT and Navy’s WindSat operate beyond their design life, many 
members of the weather and climate science communities recognize the importance of 
developing new observational technologies and strategies to meet the essential need for 
OVW information to improve hurricane intensity and location forecasts. The Hurricane 
Imaging Radiometer (HIRAD) is an innovative technology development which offers 
new and unique remotely sensed satellite observations of both extreme oceanic wind 
events and strong precipitation. It is based on the airborne Stepped Frequency 
Microwave Radiometer (SFMR), which is the only proven remote sensing technique for 
observing tropical cyclone (TC) ocean surface wind speeds and rain rates. The proposed 
HIRAD instrument advances beyond the current nadir viewing SFMR to an equivalent 
wide-swath SFMR imager using passive microwave synthetic thinned aperture 
radiometer (STAR) technology. This sensor will operate over 4-7 GHz (C-band 
frequencies) where the required TC remote sensing physics has been validated by both 
SFMR and WindSat radiometers. The instrument is described in more detail in a paper 
by Jones et al. presented to the Tropical Meteorology Special Symposium at this AMS 
Annual Meeting. 

Simulated HIRAD passes through a simulation of hurricane Frances are being developed 
to demonstrate HIRAD estimation of surface wind speed over a wide swath in the 
presence of heavy rain. These are currently being used in “quick” OSSEs (Observing 
System Simulation Experiments) with H*Wind analyses as the discriminating tool. The 
H*Wind analysis, a product of the Hurricane Research Division of NOAA’s Atlantic , 
Oceanographic and Meteorological Laboratory, brings together wind measurements from 
a variety of observation platforms into an objective analysis of the distribution of wind 
speeds in a tropical cyclone. This product is designed to improve understanding of the 
extent and strength of the wind field, and to improve the assessment of hurricane 
intensity. See httn://www.aoml.noaa.gov/hrd/data sub/wind.html . 

Observations have been simulated from both aircraft altitudes and space. The simulated 
flight patterns for the aircraft platform cases have been designed to duplicate the timing 
and flight patterns used in routine NOAA and USAF hurricane surveillance flights, and 
the spacebome case simulates a TRMM orbit and altitude. 

For the nature run from which the observations are simulated, we chose to use a 
simulation from the state-of-the-art system described by Chen et al. (2007). The storm 



simulated is Hurricane Frances (late August of 2004) using a system of nested grids with 
the innermost one having a horizontal grid spacing of 0.01 5 degrees (~1 .6 km) in 
longitude and latitude. The model is non-hydrostatic in the atmosphere with detailed 
explicit microphysics and an interactive ocean wave model (Chen et al., 2007). The 
results include a realistic eyewall, rainbands and other convective and mesoscale 
structure. (See Fig. 1.) 



Fig. 1 . Near-surface wind speed (m/s) as modeled for Hurricane Frances at 1 8z on 31 
August, 2004. Wind speed is represented by the color scale, and rain rate with solid 
contours (interval = 20 mm/hr). Only a portion of the inner-most grid is shown. 

The HIRAD sampling, as currently implemented, constructs cross-track scans with each 
gridded wind speed value so that the spatial resolution over the HIRAD field-of-view is 
of the same size as the model grid both cross-track and along track. (This resolution is of 
similar size, but not exactly the same, as the planned instrument.) Contiguous scans are 
formed along the track in a pushbroom fashion as the platform flies along. Both “Figure- 
4” and “Butterfly” patterns, which provide 2 and 3 passes through the eye, respectively, 
have been simulated for the aircraft cases and a single pass was simulated in the 
spacebome case (at 5 km resolution). Surface averaging over realistic antenna beam 
resolution cells will be added to the surface sampling simulations in the future. The 
HIRAD cross-track field-of-view is approximately ± 60 deg. which translates to a swath 
width of ~70 km at the highest aircraft altitude considered, which is 20 km for a typical 
ER-2 flight. In the satellite case the swath was approximately 2000 km for an orbit 
altitude of 450 km. In the satellite simulation chosen, the satellite gtound track passed 



through the hurricane eye, and for this example the limited modeled wind field occupied 
only ± 45 deg. of the available HIRAD field-of-view. 

Since the modeled surface wind fields are provided on 1 hour increments, time 
interpolated wind fields are used in the surface sampling. The time interpolation is done 
using storm centric 1 hour data over the duration of the aircraft flight patterns, which take 
approximately 1*4-2 hours to complete. This technique was used for all of the simulated 
observations. 

Errors that are representative of the particular instrument were added to all simulated 
observations. For HIRAD, wind speed errors were simulated using a simplified model 
that was calibrated with SFMR errors in estimating actual hurricane winds. The total 
SFMR error was separated into a surface component, rain free approximation, and an 
atmospheric component, where the standard deviation of the total error was the vector 
sum of these two. Since SFMR is nadir viewing, the standard deviation of the HIRAD 
error was modeled by applying a sec8 dependence to the atmospheric component. For 
each simulation case, this model was used in a single trial way to produce a pixel-by- 
pixel wind speed errors that were a function of modeled rain rate and viewing angle over 
the swath. Since this single trial method produces a few large errors in each simulation 
that can skew the H* Wind results, a 2 sigma limit was applied to each random error and a 
3x3 median filter was applied to the resultant wind field. 

QuikSCAT observations were simulated as follows. The objective is to simulate a 
QuikSCAT hurricane observation that is statistically representative of actual QuikSCAT 
hurricane measurements of ocean vector winds. The simulated results provide similar 
spatial sampling with reasonable error characteristics e.g., wind speed saturation for high 
wind speeds (> 30 m/s); contamination of wind speed and wind direction in the presence 
of rain, etc. The wind speed simulation incorporates a transfer function that is derived 
from the correlation of QuikSCAT L2B-12.5km retrieved wind speeds and corresponding 
QRad rain rates and H*Wind surface analysis for 10 collocated hurricane events (2003- 
2005). (See Fig. 2 for examples of the distribution of these distributions.) This transfer 
function, also called wind speed ratio (WSR), is modeled as a ratio of L2B-1 2.5km wind 
speed and H*Wind surface wind speed with respect to QRad rain rate (R), represented 
symbolically as follows: 

WSB = f(L2B 12.5, fl • Wind.R) 

Given a nature-run surface wind and rain fields, we compute the wind speed ratio; and 
the simulated QuikSCAT wind speed retrieval is the product of this WSR and the nature- 
run surface wind speed ( Wspdmuh)- 

QSCAT^ d = H'Sfl X Wsn>d M 

To simulate the QuikSCAT random wind speed retrieval error, zero-mean Gaussian noise 
is added to the simulated wind speed with the standard deviation increasing with wind 
speed. Although perfect simulation for any pixel is impossible to achieve, the simulated 



wind speed field resembles actual QuikSCAT observations both by spatial patterns and 
statistics. An example of the simulation wind speeds is shown in Figure 2 for a nature- 
run simulation of Hurricane Frances, along with actual QuikSCAT measurements for 
three hurricane over-flights. 

For clear sky and light rain, the wind direction simulation incorporates a random 
retrieved wind direction uncertainty (-18° rms). For heavier rain (> 15 mm/hr), 
QuikSCAT wind direction retrievals are cross-swath (independent of actual surface 
winds). An example of the simulated OVW for Hurricane Frances is shown in Figure 3. 






Fig.2 Simulated (upper left) and actual QuikSCAT hurricane wind speed comparisons. 







(b) Simulated QuikSCAT OVW with nature-run rain contours. 
Fig. 3 Hurricane Frances QuikSCAT simulation. 














Other wind observations simulated were: Dropsondes, GOES cloud winds, buoys, and 
ASOS. (Future observations to be added will be flight-level winds and tail-Doppler 
radar.) Dropsonde measurements were simulated from a typical flight level (3 km), and 
3-D trajectories were calculated using the model winds and assuming an empirical drop 
rate as a function of height added to the model vertical velocity. Spatially (vertically) 
correlated noise corresponding to turbulence was added to the measurements. Drops 
were made around the center of the eyewall and in the center of the eye (i.e. 6 drops per 
Figure-4 flight). See Fig. 4, which shows locations of drops and horizontal trajectories of 
the dropsondes for one flight leg. 
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Fig. 4. Nature run flight-level wind speed along with locations of sonde drops (blue 
circles) and trajectories (blue line) for one aircraft flight leg. Note that in this case, one 
of the sondes passed through the wind maximum. 

GOES cloud winds were simulated by applying a location (horizontal and vertical) 
template obtained from actual GOES winds during the time period studied, translated 
horizontally by the difference in locations of the simulated and real storms. The template 
was then used to interpolate values of the wind vectors from the nature run. The wind 
components were average vertically over 1000 meters to simulate the vertical “smearing 
effect of using clouds as tracers, and a random error of sigma 2 m/s (limited to no more 
than 20% of the wind speed value) was added to the simulated measurement. Note that 




the GOES winds were generally distant from the vortex and thus did not have a large 
influence on the H*Wind analysis. 

Buoys and ASOS measurements were simulated by interpolating the nature run winds to 
the point of measurement, and temporally correlated noise added to represent turbulence. 
However, the number and location of these measurements were such that they did not 
substantially influence the H*Wind analysis. 

H*Wind analyses were conducted for various configurations of observations. Because 
H*Wind is used by the operational community to evaluate the intensity, structure, and 
location of tropical storms, especially over ocean, the authors believe that an evaluation 
of the improvement to an H*Wind analysis by any one or set of measurements provides 
an excellent metric for assessing the relative value of various observations. 



Fig. 5. H*Wind surface wind analyses, (a — left) With full nature run information, (b 
center) with current satellites, no aircraft,; (c - right) with HIRAD satellite, no aircraft. 

Figure 5a shows the H* Wind analysis when the nature run itself is used as input, and thus 
represents the best that can possibly be done with this case with H*Wind. In comparison 
with the nature run itself (Fig. 1) the H*Wind analysis does not include much of the 
smaller-scale structure, but it does capture the maximum wind speed, the gross structure 
of the wind field including quadrants of strongest and weakest winds, and obviously the 
location of the vortex. Figure 5b shows the H* Wind analysis with no aircraft 
observations and with current satellite capabilities (QuikSCAT and GOES), and Figure 
5c shows the QuikSCAT analysis with the proposed HIRAD satellite. Since HIRAD is 
able to observe OVW through heavy rain and high winds, the vortex is clearly much 
better defined in that case. 



c. With HIRAD at 1 1 km d. With HIRAD at 20 km 


Figure 6 shows the results of H* Wind analyses with various aircraft configurations: (a) 
with SFMR; (b) with HIRAD at 3 km; (c) with HIRAD at 1 1 km; (d) with HIRAD at 20 
km. As HIRAD is flown at a higher altitude, the swath is wider and hence the advantage 
is greater. Since H*Wind does not capture high-resolution features in any case (except 
for the wind maximum itself), the lower spatial resolution due to higher altitudes does not 
degrade the H*Wind analysis. The vortex structure is much better defined for the high- 
altitude HIRAD than it is for the other cases. However, even the lower-altitude HIRAD 
provides more information than SFMR and thus results in an arguably more realistic 
analysis. For example, the maximum wind speed is closer to that of the nature run. 


